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In this work, we consider a generalised gravitational theory that contains the Einstein term, a 
scalar held and the quadratic Gauss-Bonnet term. We focus on the early-universe dynamics, and 
demonstrate that a simple choice of the coupling function between the scalar held and the Gauss- 
Bonnet term and a simplifying assumption regarding the role of the Ricci scalar can lead to new, 
analytical, elegant solutions with interesting characteristics. We hrst argue, and demonstrate in the 
context of two different models, that the presence of the Ricci scalar in the theory at early times, 
when the curvature is strong, does not ahect the actual cosmological solutions. By considering 
therefore a pure scalar-GB theory with a quadratic coupling function we derive a plethora of inter¬ 
esting, analytic solutions: for a negative coupling parameter, we obtain inhationary, de Sitter-type 
solutions or expanding solutions with a de Sitter phase in their past and a natural exit mechanism 
at later times; for a positive coupling function, we hnd instead singularity-free solutions with no 
Big-Bang singularity. We show that the aforementioned solutions arise only for this particular choice 
of coupling function, a result that may hint to some fundamental role that this coupling function 
may hold in the context of an ultimate theory. 

I. INTRODUCTION 

In the quest for the final theory, that would unify 
the gravitational interactions with their particle physics 
analogs and describe them correctly at arbitrarily large 
energy scales, the generalised gravitational theories have 
played a primary role. In their majority, these theories 
are purely geometrical in nature and involve extra dimen¬ 
sions: well-known examples are superstring theory [l|, 
the Lovelock theory 0 or the novel theories with extra 
dimensions ii- For a 4-dimensional observer, however, 
the dynamics and content of this higher-dimensional, fun¬ 
damental theory is translated into the appearance of new 
terms, describing gravity or a number of additional fields, 
in the context of the 4-dimensional effective theory. In¬ 
spired by the aforementioned theories, many variants 
of modified gravitational theories have been constructed 
over the years, and their implications for gravity and cos¬ 
mology were extensively studied. 

The most usual way to modify the gravitational in¬ 
teractions in a 4-dimensional context is via the addition 
of gravitational terms that involve higher powers of cur¬ 
vature. In the context of the heterotic superstring ef¬ 
fective theory , for instance, the Einstein term is 
supplemented by quadratic curvature terms, such as the 
Gauss-Bonnet term Rq^ or the RR term with the latter, 
however, being trivially zero for a Friedmann-Lemaitre- 
Robertson-Walker background. The Gauss-Bonnet term 
is also the first-order correction to the Einstein term 
in the Lovelock theory, that is considered as a natural 
generalisation of Einstein’s theory of gravity in a higher 
number of dimensions Q. In 4 dimensions, however, the 
Gauss-Bonnet (GB) term is a topological invariant and 
adds nothing to the held equations of the theory unless it 


is coupled to an additional held. Inspired by superstring 
theory, the usual way is to couple the GB term to a scalar 
held. By doing that, the GB term remains in the theory 
and has been shown to lead to a variet y o f new solutions: 
singularity-free cosmological solutions [3, [l^ , novel hairy 
black holes [ll|) 113 even traversable wormholes [l3| 
" solutions that were absent in the traditional General 
Relativity that contains only the Ricci scalar term (for 
a review on this type of generalised gravitational theo¬ 
ries and further references, see H). Also, since the held 
equations remain second order for the Lovelock theory in 
general and GB in particular, they are widely known to 
be ghost-free theories contrary to some brane models like 
DGP (see e.g. [l3 for a summary) or higher-derivative 
theories, which often suffer from the Ostrogradski ghost 
[13 (however, see [13 for a counterexample). These theo¬ 
ries have also been intensively studied in various realistic 
contexts: by considering solar system constraints [ll|, as 
a dark ener gy m odel with GMB and galaxy distribution 
constraints |l9j or as an inhationary model HO, [HI. 

In this work, we will focus on the dynamics of early 
universe and look for the corresponding cosmological so¬ 
lutions. During the early universe, it is natural to expect 
that a string-inspired theory would describe better its 
dynamics - that is why here we consider a 4-dimensional 
theory containing the Einstein term and a scalar held 
(one of the many of string theory) coupled non-minimally 
to gravity through a general coupling function to the 
quadratic Gauss-Bonnet term. However, we will go one 
step beyond that: we will demonstrate, in the context of 
two different models, that the presence of the Ricci scalar 
adds nothing to the dynamics of the universe at early 
times where the curvature is strong; in fact, it is the cou¬ 
pled system of scalar held and GB term that dominates 
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and governs the evolution of the universe. Motivated 
by this, we will then look for and derive, exclusively via 
analytical means, cosmological solutions with interesting 
characteristics in the context of the pure scalar-GB the¬ 
ory. For the particular case of a quadratic coupling func¬ 
tion between the scalar field and the GB term, we will 
present two classes of solutions that are distinguished by 
the sign of the coupling function constant: for a nega¬ 
tive coupling constant, the system of equations supports 
cosmological solutions that are either purely de Sitter¬ 
like or have a de Sitter inflationary phase in their past 
and a natural exit mechanism at later times; for a posi¬ 
tive coupling constant, a class of singularity-free solutions 
emerges instead. 

Although these solutions are derived for a particular 
form of the coupling function, it has been shown in the 
past that a polynomial, even coupling function, such as 
the quadratic one chosen here, shares many common 
characteristics with the actual coupling function between 
the moduli helds and the GB term in the heterotic su¬ 
perstring effective theory 0- Therefore, our present 
analysis, that proposes a new approach in deriving early- 
time cosmological solutions, will be of relevance to the 
superstring effective theory itself, but also to any other 
generalised, string-inspired theory that contains the GB 
term coupled to a scalar field. All these theories are 
bound to have in the phase-space of their solutions, the 
aforementioned classes of solutions as early-time asymp¬ 
totics. In addition, the analytical derivation of the solu¬ 
tions, in contrast to the usual numerical means employed 
systematically when the Ricci term is kept in the theory, 
allows for a more comprehensive study of the properties 
of the found solutions. The presence of the Ricci scalar 
inevitably becomes important, as the universe expands, 
and a transition between the scalar-GB-dominated phase 
and a subsequent one that prepares the universe for a 
traditional Ricci-scalar-driven cosmology will eventually 
take place. Such a transitory phase is studied at the final 
part of this work. 

Our approach is very similar to the two important con¬ 
jectures about asymptotic dynamics at early times pro¬ 
posed by Belinsky, Khalatnikov and Lifshitz (BKL). The 
two conjectures state that matter content and 

spatial derivatives are not dynamically significant 
near the initial singularity. With the same spirit, we 
will conjecture and prove in a class of models that lower- 
curvature terms in the Riemann tensor are negligible near 
the singularity compared to higher-curvature terms. 

The outline of our manuscript is as follows: in Sec¬ 
tion 2, we present the theoretical framework of our anal¬ 
ysis and the derived set of field equations. In Section 
3, through the use of a toy model, we demonstrate how 
the presence of the Ricci term in the theory merely com¬ 
plicates the analysis without changing anything in the 
actual dynamics of the solution. We reinforce this ar¬ 
gument in Section 4, where the complete theory of the 
Ricci scalar, a scalar field and the Gauss-Bonnet term 
is studied for the case of a linear coupling function. In 


Section 5, we consider the pure scalar-GB theory with 
a quadratic coupling function, and derive a variety of 
early-time cosmological solutions with attractive charac¬ 
teristics. In Section 6, we briefly consider the case of 
a general polynomial coupling function, and in Section 
7, we study an indicative transitory phase for the uni¬ 
verse as it passes from a scalar-GB-dominated phase to 
one where the Ricci scalar starts being important. We 
present our conclusions in Section 8. 


II. THE THEORETICAL FRAMEWORK 


In this work, we consider a string-inspired gravitational 
theory that contains, apart from the Ricci scalar, a scalar 
field (j) coupled non-minimally to gravity via a general 
coupling function /(</>) to the quadratic Gauss-Bonnet 
term. Such a theory is described by the following action 
functional 



(V^ 

2 


+ g /('/')-^GB 


where the Gauss-Bonnet term Rq^ is defined as 


^GB — Rp 


- ARauR^'^ + R^ . 


( 1 ) 

( 2 ) 


Since the focus of this work will be the dynamics of 
the universe at very early times where high-energy and 
strong-curvature effects are expected to be the dominant 
ones, throughout this work, we will assume that any ad¬ 
ditional distribution of matter or energy, apart from the 
scalar field, plays only a secondary role and thus will be 
ignored. 

The variation of the action CD with respect to the 
scalar field (j) and the metric tensor leads to the scalar 
and gravitational field equations, respectively; these have 
the form 

9^4>] + g f'^GB = 0 , (3) 

and 

+ , (4) 

where f' = df /dcj), and is defined as 

Pfiau^ R^ai/f5p‘^9li[0Ru]oLp‘^9a[yR!3]^PRg^]a ■ (^) 

We note that if, the Gauss-Bonnet term is altogether 
ignored in the theory, the scalar field looses its poten¬ 
tial; on the other hand, if the scalar coupling function 
is a constant, then V“^/ = 0 and the contribution of 
the Gauss-Bonnet term to the gravitational field equa¬ 
tions (CD vanishes. We may therefore conclude that (j) 
and Rqs seem to form an independent pair of quantities 
that mutually support each other. Recall that it is the 
metric that acts as potential in General relativity and it 
is coupled to the velocity, u“, in the Lagrangian, gabu‘^u^, 
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for particle motion. It should be noted that the coupling 
of the scalar field with the Gauss-Bonnet term is in the 
same general relativistic form and spirit. It is therefore 
a very desirable feature of the theory. 

We will also assume that the line-element has the 
Friedmann-Lemaitre-Robertson-Walker form 


ds^ = —dt^ + a^{t) 


dr^ 


1 — kr^ 


+ r^ {de^ + sin^ e d(p^) 


( 6 ) 

that describes a homogeneous and isotropic universe with 
a scale factor a{t) and spatial curvature k = 0, ±1. For 
the above metric ansatz, the Gauss-Bonnet term takes 
the explicit form 


R%s = 24{H^ + ^)iH^ + H), (7) 


where H = d/a is the Hubble parameter and the dot 
denotes derivative with respect to time. Using the above 
geometrical quantities, the scalar m and gravitational 
field equations (jl]) reduce to the following system of three, 
ordinary but coupled, differential equations 


^ + 3H^-3f{H^ + ^)iH^+H)=0, ( 8 ) 

3il + Hf){H^ + ^) = ^, (9) 

2(1 -b Hf)iH^ + H) + {1 + f){H^ + A) = -^ . (10) 

In the following sections, we will look for cosmological 
solutions supported by the above set of equations both 
in the presence and in the absence of the Ricci scalar. 
We will start with a simple toy model, that will demon¬ 
strate the role of the Ricci scalar in the context of the 
complete theory; the same task will then be performed 
in the framework of a less restricted model. Using the 
derived insight of when the Ricci scalar may be ignored 
from the theory, we will then search for solutions where 
the Ricci scalar is negligible compared to or of the same 
order as the Gauss-Bonnet term. 


III. A TOY MODEL 


while for fc = — I we obtain a linearly expanding universe 
with a(t) = At + B and a singularity at a finite time. 

We will thus search for solutions with 0^0. Since 
Rqb = 0, the scalar equation is easily integrated once to 
give 

= ( 12 ) 

where C an arbitrary integration constant. The solution 
for the scale factor is most easily given by the constraint 
Rqb — 0 itself, or equivalently 



where the two multiplying factors could be independently 
zero or not. In fact, it is only for the second choice, d = 0, 
that a non-trivial solution for the scalar field is allowed. 
Then, we easily find that a{t) = At + B again, but now 
this solution holds for all values of k. Taking the sum of 
Eqs. ds-cni), we find the constraint 

/ + 3/-+4 = 0, (14) 

a 


with solution 


fit) = - 


Cl 


2A{At + H)2 


Bt 

~A 


A 

2 


C2, 


(15) 


where Ci ^2 are again integration constants. The scalar 
field itself can easily be found from Eq. (IT^ to have the 
form 


cj){t) = 


C 

2A{At + B)'^ 


+ D. 


(16) 


From a field-theory point-of-view, it would be much 
preferable to express the coupling function / in terms 
of the field (j) instead of the time coordinate. Comparing 
the above two equations, and upon a convenient choice 
of the integration parameter D, i.e. D = 0, we see that 
we can write 


fW = fl</ + J + f3, 


(17) 


In this section, we will look for solutions of the set 
of Eqs. (I51)- (ITU1) on which we impose the constraint of 
the vanishing of the Gauss-Bonnet term, Rqb — 0- The 
(j) = const, is then a solution of the scalar field equation 
which is trivially satisfied. But, the Einstein’s equations, 
since now (/> = (/ = Q, demand that the following two 
constraints 


where fi are constants given in terms of (A, B, C, ci, C 2 ). 

Let us now ignore the presence of the Ricci scalar in 
the theory. Since the scalar-field equation dH) and the 
constraint (US remain unaltered, both the solution for 
the scalar field (fTBl) and the linearly-expanding form of 
the scale factor are still valid. However, the Einstein field 
equations now take the simplified from 


A!_A=o, and d = 0 (II) 

should be simultaneously satisfied. As it was shown in 
0 , for k = 0,4-1, the above equations lead to a static 
universe with an arbitrary or infinite radius, respectively. 


3{k + d'^) -a 
a? a 



( 18 ) 




a a 


(19) 
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which, when combined with each other, lead to the con¬ 
straint 


/ + 3/-=0, (20) 

a 

for the coupling function. The above can be easily inte¬ 
grated and expressed in terms of the scalar field to obtain 

/(<^) = /o<^ + /3. (21) 

Therefore, assuming that the Ricci scalar can be ignored, 
we obtain the same solutions for the scalar field and scale 
factor as in the presence of it, with the only change ap¬ 
pearing in the expression of the coupling function that 
now takes a simpler form. Looking more carefully, the 
two expressions dnl) and (I20|l are equivalent in the limit 
where the scalar field takes very large values; from Eq. 
(nsj, this happens as we approach the initial singularity, 
{At + B) ^ 0. Therefore, as expected, in the presence of 
the quadratic Gauss-Bonnet term in the theory the Ricci 
scalar adds nothing to the dynamics of the universe at 
the very early-time limit. 


IV. THE COMPLETE THEORY WITH A 
LINEAR COUPLING 


In this section, we will attempt to reinforce the con¬ 
clusion drawn in the previous section regarding the role 
of the Ricci scalar in the early-universe cosmology but in 
the context of a more realistic set-up. We will therefore 
look for physically-interesting solutions following from 
the complete set of Eqs. (l8])- (fT0l) under the assumptions 
that the Gauss-Bonnet term is not zero and that the func¬ 
tion f{4>) is a non-trivial function of the field tp- Then, 
using the relations 

/ = />, / = /'>2^/>, (22) 


and adding Eqs. m and (iinD, we end up with the con¬ 
straint 


(i/2 + ^) [4 + f{4, + 3i/<^) + /'V2 


+ 2{H^ + H){l + fH,j))=0. (23) 


We may now use Eqs. ([8]) and Q to replace the combi¬ 
nation {(j) + 3H(p) and (j)'^ in the second and third term, 
respectively, of the above equation. Then, we arrive at 


ignoring the Gauss-Bonnet term from the theory. The 
particular case of n = 1 will be studied in this section 
while the case with n = 2 will be considered in Section 

5. 

Let us therefore focus on the scalar equation (|8l) and 
use that f{cj)) = Xcj); then, it can be brought to the form 


= 3Ad(fc + a2), (25) 

at 

which can be straightforwardly integrated to yield the 
relation 


C Xa {3k + a^) 


(26) 


with C an integration constant. The second term in the 
above expression, proportional to the coupling constant 
A, is clearly the one associated to the Gauss-Bonnet term, 
while the first one appeared also in Eq. m in the ab¬ 
sence of a potential. As we are now interested in studying 
the properties of solutions arising in the presence of the 
quadratic GB term, we will set for simplicity (7 = 0, and 
make a comment on the role of a non-vanishing value of 
C later in this section. 

Eor the case of linear coupling, Eq. (l24l) is also sim¬ 
plified since /" = 0. Moreover, we may use Eq. ([5]) to 
replace the combination {l + f'H(j>) in terms of . Then, 
Eq. takes the simpler, more explicit form 


2^2 


3{k + d^) 


4 + A2^(fc + a2) 


.3 i2 


= 0 . 


(27) 


Now, keeping only the second term in Eq. (E5)) and using 
this relation to substitute cj) in the equation above, the 
latter constraint is finally rewritten as 


3{k + a2)2 


4 + A2^(yfc + a2) 


+ A2^(3fc + a2)2 = 0 . 


(28) 

The above equation does not involve </> anymore, and it 
can be integrated once to give 


3a^ = -A2 


/ 13ka? 


50^ 2k^ \ 

~ ^ k + a? ) 


+ Cl , 


(29) 


with Cl an integration constant again. Unfortunately, 
for k ^ Q, the above cannot be easily solved for a to 
yield, via another integration, the form of the scale factor. 
However, for the case of a flat universe {k = 0), we easily 
write that 


4 ( 7/2 +A) + 2(i/2 + i/) + 6/"(/i2 + ^y 


(1 + f'H^) 



(30) 


+ 3/'2(i/2 +A^yi/2 + ii) = 0. (24) 

In this work, we will consider a polynomial form for 
the scalar coupling function, i.e. /(</)) = A^", where A a 
constant and n an integer. The case with n = 0 results 
into a constant coupling function and is equivalent to 


with solution 

= + (31) 

In the above, F{a,b^c;x) stands for the hypergeometric 
function whose convergence demands that a{t) < a^ax = 


a{t) F 


I I 5 3a^{t) 
4’ 4’ 4’ 


Cl 
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FIG. 1: The scale factor a{t) in terms of time, for a lin¬ 
ear coupling function and A = 0.1, and for the values ci = 
0.1, 0.5, 1, 2, 5 (from bottom to top). 


(ci/3)^/"^ - from Eq. (j29l) and for fc = 0, it follows that 
Cl must be indeed a positive constant for a real solution 
for the scale factor a{t) to exist. As a result, the scale 
factor is bounded from above with the constant ci deter¬ 
mining its upper value. On the other hand, in the limit 
a —>■ 0, the hypergeometric function goes to unity, and 
the relation between the scale factor and the time coordi¬ 
nate becomes linear; this of course signals the existence 
of a singularity at a finite value of the time-coordinate. 
Figure [T] depicts the proHle of the scale factor as it starts 
from the initial singularity, follows an increasing phase 
and eventually reaches its maximum value determined 
by Cl. The value of the coupling parameter A affects the 
slope of the uprising curve at early times. 

Let us make, at this point, a comment on the role of a 
non-vanishing value of the parameter C. If we keep both 
terms in Eq. (I2S1) and follow a similar analysis as above, 
then, for fc = 0 again, we obtain, instead of Eq. (l30l) . the 
following equation 

C"^ t; \2 

3a^ = — - 2XCa - — -b ci . (32) 

20 '^ 2 

The above is again difficult to solve for d and integrate 
further, however, we may study particular limits. For a 
non-vanishing value of C, ignoring the first two terms, as 
we did in the analysis above, would still be justified in a 
regime where d takes very large values. Drawing expe¬ 
rience from the known cosmological solutions, we know 
that the rate of expansion of the universe is the largest 
close to the initial singularity; as we are indeed interested 
in the early-universe dynamics, the role of C is therefore 
negligible and our assumption was justified.^ 


^ For completeness, let us note that in a regime where a would 
be very small and thus the first term on the right-hand-side of 
Eq. II32II would dominate instead, the corresponding solution for 
the scale factor is given by the form a(t) ~ (At + . For 

intermediate values of a, where the second term in Eq. (|32} 
would dominate, the scale factor behaves as a(t) ~ (t -|- . 


Looking for further support to our argument that the 
presence of the Ricci scalar adds very little to the dynam¬ 
ics of the early-time cosmological solutions, we will now 
adopt a different approach: we will focus on the early- 
time limit, and ignore from the beginning the Ricci scalar 
from the field equations. In the absence of R, the grav¬ 
itational equation (TTSl) can be straightforwardly solved 
for (/) to yield 

<^ = 6/'ij(ifc2 +A) , (33) 

even for a general coupling function f{4>) and without any 
need for integration. If we then follow a similar analysis 
as at the beginning of this section, i.e. take the sum of 
Eqs. m and m, and use the scalar equation ([5]) and 
Eq. (1331) in there, we obtain the new constraint 



0. (34) 


The quantity f is not allowed to be zero since then the 
Gauss-Bonnet term would be eliminated from the theory. 
The same holds for the combination {H^ + k/a^) since 
this term is part of both the Gauss-Bonnet expression 
and (p, Eqs. © and (1551) respectively. Therefore, it is the 
expression inside the square brackets that should vanish 
instead; for a general polynomial form of the coupling 
fucntion, /{(j)) = this may be re-written as 


12a 


(fc -b 5a^) d -\ -^ (fc -b d^y \n(n — 1) (j) 


xn —2 


0. (35) 


Therefore, for the linear case with n = 1, the second 
term in the above equation vanishes which leaves us with 
a much simpler constraint: its general solution is a linear 
function of time, i.e. a{t) = At+B. This is in accordance 
to the behaviour obtained in the first part of this section 
when the early-time limit of the complete solution (1511) 
was considered. The dashed lines in Fig. [T] represent the 
linear-in-time solution for the scale factor, derived above 
in the absence of the Ricci scalar in the theory - the 
agreement with the early-time behaviour of the complete 
solution is more than evident. Therefore, we conclude 
again that, by ignoring the presence of the Ricci scalar, 
we significantly simplify the analysis and still obtain ex¬ 
actly the same early-time cosmological solution. 

As we anticipate from the above, the solution for the 
scalar field may be derived, in the same early-time limit, 
from either Eq. (EH) with C = 0 or from its simpler 
analog (1551) : in both cases, we easily obtain that at early 
times 


= (f)o 


3A (fc-b A2) 
{At -b By 


(36) 


where A, B, and (j)o are arbitrary integration constants. 
The above describes a decaying scalar field as the uni¬ 
verse expands with a singular behaviour at the initial 
singularity. 
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V. THE SCALAR-GB THEORY WITH 
QUADRATIC COUPLING 


In this section, we address the case of the quadratic 
coupling function of the scalar field to the Gauss-Bonnet 
term, / = X(j>^. The constraint (l24l) . derived in the con¬ 
text of the complete Einstein-scalar-GB theory, is valid 
for an arbitrary /(^i), therefore it could be applied for 
this case, too. However, the scalar equation (I5|) now can¬ 
not be easily integrated, and as a result the set of field 
equations cannot be decoupled. Despite our persistent 
efforts, no way forward could be found via analytical cal¬ 
culations. 

Nevertheless, if one is interested strictly in the early- 
universe dynamics, the results of Sections 3 and 4 point 
towards ignoring the Ricci term from the very beginning 
in order to simplify the analysis and increase the chances 
of deriving viable cosmological solutions via analytical 
means. Equation (I35L valid in the context of the pure 
scalar-Gauss-Bonnet theory, has been derived for a gen¬ 
eral polynomial coupling function and it can be applied 
directly for the quadratic case with n = 2. In that case, 
the 0-dependence disappears and the constraint becomes 
a differential equation only for the scale factor a{t). It 
can be conveniently separated as follows 

(fc-|-5d^)d 24Ad 

d(/c-hd2)2 = ^ ^ 

leading eventually, after integrating both sides with re¬ 
spect to time, to the relation 


12A 



2 

{k + a?) 


+ Ci. 


(38) 


In the above, Ci is an arbitrary integration constant. For 
k ^ 0, the above relation is of a transcendental form and 
thus impossible to solve for d. We are thus forced to 
consider again the case of a flat universe, in which case 
we obtain the simple differential equation 


2 d2 ^ a2 


(39) 


In fact, a variety of cosmological solutions with interest¬ 
ing characteristics may be derived from the above simple 
equation depending on the values of the integration con¬ 
stant Cl and the Gauss-Bonnet coupling parameter A. 
Below, we present a comprehensive analysis of all the 
ensuing solutions. 


A. The case with Ci = 0 


If Cl = 0, then from Eq. (1391) it is clear that solutions 
exist only for A < 0. In this case, a simple integration 
yields the solution: 


a(t) = ao exp 



(40) 




FIG. 2: The scale factor a(t) and scalar field 0(t) ver¬ 
sus time for the cosmological solutions with Ci = 0 and 
|A| = 0.1,0.2,0.5,1,2,10. For simplicity, the initial values 
ao and 0o have been normalised to unity. 


The theory admits both increasing and decreasing solu¬ 
tions for the scale factor with respect to time - choosing 
the positive sign, we obtain an exponentially expanding 
universe with no singularities at finite values of the time 
coordinate. Note that no self-coupling potential V(^) 
needed to be introduced for the scalar field or tailored in 
any ad hoc way in order to obtain inflation. Instead, 
it is the coupling of the scalar held to the quadratic 
Gauss-Bonnet term that provides a potential in the most 
natural way and supports inhationary solutions in the 
early universe. Such couplings arise naturally both in 
the context of string-inspired or gravitationally modihed 
theories, and we anticipate that they should all contain 
similar inhationary solutions - by insisting, however, in 
keeping the Ricci scalar in the analysis, these solutions 
were missed. 

The form of the scalar held, in turn, may be easily 
found via Eq. (1551) . that now takes the form 



(41) 


Employing the solution for the scale factor (HOI) found 
above (with the positive sign), the scalar held is found to 
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have the form 


The scalar field is also everywhere well-defined and is 
decaying exponentially from an arbitrary initial value 
to zero. The profiles of both the scale factor and the 
scalar field are shown in Figs. [Ha,b), respectively, for the 
values |A| = 0.1, 0.2, 0.5,1, 2 and 10; we observe that the 
smaller the value of the coupling constant A, the faster 
both quantities evolve with time. 

The effective potential of the scalar field receives 
contributions from both the coupling function and 
the Gauss-Bonnet term, and has the form 14// = 
—/((/))i?Qg/8 ~ cjP'/\X\. Further aspects of this de 
Sitter solution were studied in a previous short work 
[^ . In there, we showed that the necessary number 
of e-foldings follow easily without the need for assuming 
trans-planckian initial values (j)Q for the scalar field. As a 
result, the effective potential, although of a similar form 
to that of chaotic inflation [^, remains always bounded. 
Its dependence on the coupling constant A allows also 
for its value, at the time of inflation, to be large enough 
so that it dominates indeed over the other constituents 
of the universe. Finally, its quadratic form places it in 
the group of inflationary models that are still compatible 
with the current observational constraints [Slj. 


B. The case with Gi > 0 and A < 0 

We now assume that the coupling parameter A is again 
negative as before, but we allow the integration constant 
Cl to take positive values. This class of solutions first 
appeared in [2^ but we include it also in the analysis 
here for completeness and for providing additional math¬ 
ematical details left out in the previous work. In this 
case, Eq. (1551) can be rewritten as: 

/ + = + (43) 

where = I2|A|/Ci. By setting a = v tanw and using 
standard techniques of integration, the above equation 
eventually leads to the result 

/—;: - - ( \/a?‘ + v'^ — v\ I 5 , , 

^-j + 

(44) 

For z/ = 0, i.e. in the absence of the Gauss-Bonnet 
term, we smoothly recover the linearly expanding solu¬ 
tion found in Section 3. For zz 0, Eq. (l44l) is of a 
transcendental form and thus cannot be solved for a as 
a function of t. Its behaviour in terms of time is nev¬ 
ertheless depicted in Fig. [3J choosing the (-l-)-sign, we 
find that the scale factor is expanding with time, first 
with a much faster pace and later with a significantly 
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FIG. 3: Cosmological solutions with Ci > 0 and A < 0, and 
for 1 / = 0.5, 0.8,1,2,3. 


slower one. As z/ increases, the point where the scale fac¬ 
tor vanishes, i.e. the initial singularity, moves gradually 
towards larger negative values of the time coordinate; on 
the other hand, for larger values of time, the variation of 
zz affects much less the rate of expansion. We may derive 
analytic asymptotic solutions describing the behaviour of 
the scale factor at early times and later times, by taking 
appropriate limits of the complete solution (|44ll . To this 
end, we first consider the limit a 0, in which case we 
recover the pure de Sitter solution (1401) found in the pre¬ 
vious subsection; we thus conclude that the families of 
solutions with A < 0 and with Ci = 0 and Ci > 0 are 
smoothly connected in the phase-space of the theory. On 
the other hand, expanding for large values of a, i.e. for 

^ z/^, we obtain a linear function of time for the scale 
factor. Therefore, in this class of solutions, the universe 
interpolates between a pure de Sitter solution at very 
early times, triggered by the Gauss-Bonnet term, and a 
linearly-expanding Milne-type phase at later times when 
the effect of the Gauss-Bonnet term starts to wear off. 
In this way, these solutions may readily accommodate an 
early inflationary phase with a natural exit mechanism 
at later times. 

The solution for the scalar field should in principle fol¬ 
low from Eq. (l33l) . however the absence of an explicit 
solution for the scale factor in terms of time complicates 
the necessary integration. However, we may find an im¬ 
plicit expression for the scalar field in the following way: 
if we use Eq. (1331) . then Eq. (1351) may be re-written, for 
a general polynomial coupling function, as 


- 2(n 



d (fc + 5d^) 
d{k + a'^) 


(45) 


which, upon integration of both sides with respect to 
time, gives the constraint 




Co 

d{k + d^)2 ’ 


(46) 


Specialising for the case k = 0 and n = 2, and using Eq. 
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(l39ll to replace a, we finally obtain 




(47) 


Using the asymptotic behaviour of the scale factor at 
early times, i.e. the pure de Sitter solution (|40)) . we may 
easily conclude that the scalar field also assumes its expo¬ 
nentially decaying form of Eq. (H^ . On the other hand, 
for ^ z/^, the scalar field reduces to a constant. 


C. The case with Ci < 0 and A < 0 



-8 -6 -4-2 0 2 

t 


We will next assume that A < 0 again, but that Ci < 0 
now. From Eq. dsni), we may see that any solutions that 
would follow in this case will not allow the scale factor to 
grow indefinitely but only up to a maximum value, oth¬ 
erwise the rate of expansion would become imaginary. 
Equation (I39L for this choice of parameters, can be al¬ 
ternatively written as: 


FIG. 4: Cosmological solutions with Ci < 0 and A < 0, and 
for = 1, 2, 3, 4, 5 from left to right. 

zero. As in the previous case, the scalar field remains 
finite and diverges only at the initial singularity where 
the scale factor vanishes. 



{t + to), 


(48) 


D. The case with Ci > 0 and A > 0 


where = 12|A|/|C'i|. We now set a = f' sinw, and by 
integrating once we find the result 

(49) 

The profile of the scale factor in terms of time, as this 
follows from the above relation for the positive sign, is 
depicted in Fig. 01 The early-time behaviour of the above 
solution is very similar to the one found in the previous 
subsection: a singularity, where a = 0, is again present 
but this is reached at increasingly large negative values of 
the time coordinate as {> increases; the expansion of the 
relation (|4^ in the limit a —>■ 0 leads once again to the 
pure de Sitter solution (HOI) and to an inflationary phase 
for the universe. However, the quantity v is now the 
maximum allowed value of a{t), and the universe stops 
increasing after this point. We thus conclude that, in the 
phase space of the solutions of the theory, the pure de 
Sitter solutions for A < 0 and Ci = 0 branch off to two 
families of cosmological solutions with totally different 
behaviour at larger times: one that allows for indefinite 
expansion of the universe for Ci > 0 and one where the 
scale factor reaches a ceiling for Ci < 0. 

The solution for the scalar field follows again from Eq. 
(H51) and in this case is given by 


,2 = Co 


(50) 


Again, in the early-time limit where a —> 0, the scalar 
field adopts the exponentially decaying form of Eq. . 
On the other hand, for d} —>■ z>^, the scalar field goes to 


We now reverse the sign of the coupling parameter A 
and allow it to take only positive values. In that case, 
the choices Ci < 0 are not allowed, and mathematically 
consistent solutions may be derived only for Ci > 0. In 
this case, Eq. dSi) can be rewritten as: 

/ ^Va2-z>2 = ±y^(t + to), (51) 

where iP' = YlXjCx. We now set a = v/ cosm, and upon 
integrating, we find the relation 

\/a2 - - z^arccos ^ 2^ 

For z/ = 0,i.e. A = 0, we may easily see that we go back to 
a linearly expanding solution with an initial singularity 
at a finite value of the time coordinate. However, for 
z> ^ 0, i.e. in the presence of the Gauss-Bonnet term, 
the situation is radically different: the square-root on 
the left-hand-side demands that > z>^, therefore v 
becomes the smallest allowed value of the scale factor in 
this model, and no singularities are allowed to arise. The 
same conclusion follows from the inverse cosine function 
whose argument should also satisfy the inequality — 1 < 
z//a<l, ori/<ain our case since both v and a are 
positive-definite. If we expand the solution (15^ for values 
of a close to its minimum value v, we find the asymptotic 
solution 


ait) ~ [I + (At + H)2/3] , (53) 

which shows its regular behaviour for any finite values of 
the time-coordinate. In Fig. O we plot the behaviour of 
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FIG. 5: Cosmological solutions with Ci > 0 and A > 0, and 
for V = 0.1, 0.5,1,1.5, 2 from left to right. 


In the present work, we have indeed managed to de¬ 
rive an analytical singularity-free solution for the even 
coupling function /(0) = in total agreement with 
the aforementioned argument. Note, however, that the 
analysis of 0 took into account the presence of the Ricci 
scalar in the theory thus forcing the authors to perform 
numerical integration in order to find the sought solu¬ 
tions. Here, by ignoring the Ricci scalar, we have man¬ 
aged to demonstrate in a very simple way the emergence 
of singularity-free solutions and to find their exact ana¬ 
lytical form. We believe that this result solidifies even 
more our argument that the Ricci scalar may indeed be 
ignored in the early-time limit with no effect in the dy¬ 
namics of the universe and that, in fact, one should do 
so, in order to derive elegant, analytical solutions. 


the scale factor as a function of time: the different curves 
correspond to different values ofz/ = 0.1,0.5,l,1.5,2 from 
left to right. We observe that as v increases the curve is 
moving towards larger values of a and thus away from 
the singularities. For large values of a, the dependence 
becomes linear again in terms of time - note that this 
behaviour is common for the two solutions arising for 
Cl > 0. 

The implicit solution for the scalar field in this case is 
given by 


<^" = Co 


/2Ciy/' (a2-p2)5/2 
V 5 j 


(54) 


Since the minimum value of the scale factor is v, the 
above expression for the scalar field remains always finite. 
It starts from a zero value at very early time, it follows an 
increasing profile reaching eventually a constant value. 

Let us comment at this point on the implications of the 
existence of this particular class of solutions. Singularity- 
free solutions arising in the context of string-inspired the¬ 
ories, or modified gravitational theories in general, are 
always important as they provide an alternative to the 
cosmological solutions of the classical theory of General 
Relativity that always possess an initial (or final) singu¬ 
larity. Solutions free from the initial singularity have 
been derived, for instance, in the context of the het¬ 
erotic superstring effective theory In there, the cou¬ 
pling function between the scalar field and the Gauss- 
Bonnet term is given by the highly non-trivial expres¬ 
sion /(())) = \n[2e^ri(i^)], where ry is the Dedekind func¬ 
tion. However, in |10l | analytical arguments were pre¬ 
sented that showed that singularity-free solutions may 
arise in the context of a similar theory where the coupling 
function has a much simpler form, namely /(</>) = A^", 
where n was strictly an even, positive number - specific 
solutions were found numerically in the context of the 
same work. The reason was that these two different forms 
share in fact three important characteristics: they are 
invariant under the change ()>—>■—(/), they have a global 
minimum and asymptotically they tend to inhnity. 


VI. THE CASE OF A GENERAL POLYNOMIAL 
COUPLING FUNCTION 

In this section, we will briefly address the case of a 
general polynomial coupling function, f{(j)) = Xcj)^. Both 
Eps. (I35|l and (gSl) have been derived for this case, and 
thus can be straightforwardly used for our purpose. Solv¬ 
ing Eq. (I46p for the scalar field (j) and substituting in Eq. 
(I35L we obtain the constraint 

(A:-h5d2)dd^^a3-f 12r(/c-kd2)"/("-i) =0, (55) 


with 

r = An(n-l)C'J^ . (56) 

Specializing to the case of a flat universe with A: = 0, 
we obtain 

5dd^lT^a^-k 12r = 0. (57) 


We can easily check that the above equation admits de 
Sitter-type solutions only for particular values of the in¬ 
teger n. Setting a{t) = exp(iJoA), the aforementioned 
equation becomes 


exp 


'5 (n- 2) 
2 {n- 1) 



-k i 2 r = 0 . 


(58) 


According to the above, consistent de Sitter-type solu¬ 
tions arise only for n = 2, and for T < 0 or equivalently 
A < 0, in accordance to the results of Section 5. 

If we integrate Eq. dSll) with respect to time once, we 
find the relation 

d^> = 4+^, (59) 

where we have defined the parameters 


_ 3(n ~ 6) n — 6 

5(n — 1) ’ ^ 2(n — 1) ’ 


(60) 


and P is an arbitrary integration parameter. We observe 
that for large values of the scale factor, the term a/a? 












10 



t 

FIG. 6: Early-time cosmological solutions in the scalar-GB 
theory with a general polynomial coupling function = 

and for various values of n. 


on the right-hand-side of Eq. (l55|) can be ignored. Then, 
for /3 = 0, the universe reaches asymptotically a static 
Einstein-type solution, while for /3 0, we obtain a uni¬ 

versal behaviour of the form a{t) = At+ B independently 
of the exact value of n. 

On the other hand, in the early universe where a is 
small, we can approximate Eq. (1591) by the following 
equation 


a 


(61) 


that upon integration, and for a > 0 or equivalently A > 
0, gives 


n-6 

a{t) = (^At + ^ ^ 02 ) 

Contrary to the universal behaviour characterising the 
large-a limit found above, the dynamics of the universe 
in the small-a limit, i.e. in the strong-curvature regime, 
strongly depends on the exact form of the coupling func¬ 
tion and thus on the value of the integer n. Note that 
due to the expression of the solution (15^ . the case with 
n = 2 - that was studied in detail in Section 5 - should be 
excluded; therefore, the solution (15^ describes the early- 
universe dynamics for solutions obtained in the context 
of the scalar-GB theory with a polynomial coupling func¬ 
tion with n ^ 2. For instance, for n = 1, we correctly 
recover the linear early-time solution a{t) = At + B ob¬ 
tained in Section 4. In general, all cases with n < 2 or 
n > 6 lead to expanding cosmological solutions with an 
initial singularity at a finite value of the time-coordinate, 
whereas all cases with 2 < n < 6 describe contracting 
cosmological solutions with a final singularity at asymp¬ 
totically infinite time. Indicative cases of the above so¬ 
lutions are depicted in Fig. |6l Clearly, it is only for 
the case with n = 2 and A > 0 - depicted also in Fig. 
m for completeness - that an expanding singularity-free 
solution arises in the theory. 


VII. A TRANSITORY SOLUTION 


Looking at the gravitational equations dU-dlol) , we 
conclude that, by ignoring the Ricci scalar in the previous 
sections, what we have actually assumed is the following 
conditions to hold 

/^»1, /»1, (63) 

even for a general coupling function /(</>). In sections 3 
and 4, we demonstrated, for two different choices of the 
coupling function, that the Ricci scalar can indeed be 
ignored in the early-time limit and therefore the above 
conditions hold. 

In order to study how the system will transit from 
the very early-time epoch, where the Gauss-Bonnet term 
dominates, to a subsequent era, where the Ricci scalar 
starts becoming important, here we will assume that the 
following constraint holds between / and the scale factor 


/ - = Co . (64) 

a 

In the above, cq is an arbitrary constant. We will thus 
assume that there is an era where the Gauss-Bonnet 
term starts giving a contribution to the gravitational held 
equations that is of the same order as that of the Ricci 
scalar. We will hrst consider the case where Cg is still 
much larger than unity, therefore, Ricci scalar can still 
be ignored. Then, we will assume that cq ~ 1, in which 
case both gravitational terms should be taken into ac¬ 
count in the held equations. Finally, the limit cq <C 1 
will be considered, valid in a era where the Ricci scalar 
will be the dominant term. Although the exact study of 
the evolution of the system demands a numerical anal¬ 
ysis, we believe that the following study will give us a 
feeling of whether such a transition is possible. 

Assuming that the constraint (IM)) holds and that cg ^ 
1, the Ricci scalar may be ignored and the corresponding 
held equation (fTSl) is re-written as 




+ a^) 


(65) 


Also, differentiating Eq. (IMl) once with respect to time, 
we obtain 


/ = /(---). ( 66 ) 

\a a J 

When the above and Eq. (1551) are used into the second 
gravitational equation m, the latter also takes the new 
form 

^ ^ a (fc- d^) 

a a {k + d^) 

The above can easily be integrated once, using elemen¬ 
tary methods, and leads to the relation 

(fc -I- d^) a'^ = C d. 


( 68 ) 









11 


Rearranging the above and integrating once more, we 
obtain the integral equation 


2a^da 


C + eVC'2 - 4fca8 


— t tQ . 


(69) 


where e = ±1. For k ^ 0, this arbitrary sign denotes the 
presence of two distinct branch solutions that have the 
form 


8ea^ 5 13 4/ca®^ , eVC^ - 4fca8 - C 

^ 8’ T’ ~C^) ^ 6^3 - ^ + ^0 , 

(70) 

where F{a, b, c; x) is again the hypergeometric function. 
Above, we have assumed that the integration constant 
C is positive - in fact, there is a symmetry between the 
signs of C and e, therefore we are allowed to fix the sign 
of one of the two parameters. Depending on the values 
of fc = ± and e = ±1, the above expression describes 
a variety of smooth cosmological solutions with various 
asymptotic behaviours as a —>■ 0 or a —>■ oo. For instance, 
for e = — 1 and for both values of fc = ±1, solutions arise 
that do not possess any singularity at a finite value of 
the time-coordinate. 

However, here we are mainly interested in studying the 
transition of the system through the different epochs, and 
for this reason we may simplify our analysis by setting 
fc = 0. Then, from Eq. (I69L we must necessarily have 
e = 1, and a simple integration leads to the power-law 
solution 



FIG. 7: An indicative transitory solution with the value of 
the parameter co parametrising the weight of the GB term 
in the theory - the curve in the middle corresponds to the 
limiting case with co — 1. 


Then, using both (175]) and (1661) in Eq. m, we take the 
alternative constraint 

. (fc-|-d^) , 2d d (fc-hd^) „ 

4(l-|-co)- 2 -l"(l+co)-Co- 72 -— 0, (76) 

a a 

or, after a little bit of algebra, 


4(1+ co)- 
a 


a [cofc — (2 -I- co)d^] 
d (fc -I- a?) 


(77) 


a{t) = {At + Bf/^, (71) 

where A and B are again integration constants. The 
solution for the scalar field may easily follow from Eq. 
(1551) upon using the above solution for the scale factor; 
it has the form 

= (/)o ± \n{At B). (72) 

5 

For these simple forms of the scale factor and scalar field, 
one may determine the form of the coupling function 
integrating Eq. (l64l) . we obtain 

/(^) =/o + (At -I- H)^ , (73) 

which may be alternatively written as 

/(<(>) = fo + fi exp • (74) 

Therefore, the above special solution actually corre¬ 
sponds to a particular choice of an exponential coupling 
function. 

We now assume that, as the time goes by, the effect of 
the Gauss-Bonnet term begins to diminish, and the point 
is reached where cq — 1. Then, the contribution of the 
Ricci scalar must be restored, and in that case Eq. (0 
leads to the result 

0^ = 6(l+e,)lD^. (75) 


The above can be integrated once with respect to time 
to give the relation 

(fc-f d2)a4 = C'd"°Ai+=o)^ (78) 

where C an integration constant. Once again, for fc ^ 
0, the analytical calculation is extremely difficult as 
the above equation exhibits a non-algebraic dependence. 
Therefore, we set fc = 0 and, upon integrating once more, 
we find the solution 

a{t) = {At + B)(2+“o)/(6-e5co) ^ 

From Eq. GSD, we see that, for fc = 0, 1 -I- cq > 0, 
therefore the power of {At -|- R) in the above expres¬ 
sion is also positive. Therefore, the above describes an 
expanding universe with an initial singularity emerging 
when At + B ^ 0 i.e. at a finite time-coordinate. For 
Co ^ 1, the above power reduces to 1/5, as expected, in 
accordance with Eq. (ED. In the limit cg 1, the power 
becomes 1/3. For cq — 1, the power interpolates between 
1/5 and 1/3. In Fig. [71 we plot the two indicative cases 
with Co 1 and cg <C 1 as well as the limiting case with 
Co — 1 that lies in between. Although the initial sin¬ 
gularity in this particular solution cannot be altogether 
avoided, the presence of the Gauss-Bonnet term works 
towards making the singularity softer. As the universe 
expands, the rate of expansion increases due to the de¬ 
crease of the effect of the Gauss-Bonnet term reaching 
eventually its highest value, taken in the context of the 
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pure Einstein-scalar field theory, when the Gauss-Bonnet 
becomes negligible. Of course, a complete cosmological 
analysis would demand the introduction of additional in¬ 
gredients in the universe, such as the radiation or matter 
energy density, and the detailed study of the sequence 
of different eras - however, this is beyond the scope of 
the present work that investigates the role of the Gauss- 
Bonnet term in the very-early-universe dynamics. 


VIII. CONCLUSIONS 

In this work, we have considered a generalised gravi¬ 
tational theory that contained, apart from the Einstein 
term, a scalar field and a higher-curvature - the quadratic 
- Gauss-Bonnet term. Both of these additions are met 
in the superstring effective theory as well as in a variety 
of modified gravitational theories considered in the liter¬ 
ature over the years. Theories of this type, with either 
stringy-like or more general coupling functions between 
the scalar field and the Gauss-Bonnet term - a necessary 
feature in order for the GB term to remain in the the¬ 
ory - have shown to lead to novel gravitational solutions. 
Here, we have focused on the early-universe dynamics 
and demonstrated that a simple choice of the coupling 
function and a simplifying assumption regarding the role 
of the Ricci scalar can lead to new, analytical, elegant 
solutions with interesting characteristics. 

Starting from the latter basic element of our analysis ~ 
the role of the Ricci scalar in the early-universe dynamics, 
a toy model, that was considered in Section 3, hinted to 
the fact that the presence of the Ricci scalar in the theory 
merely makes the analysis much more complex without 
affecting the actual solutions for the scale factor of the 
universe and the scalar field. This hint was changed to a 
certainty when, in Section 4, the Einstein-scalar-GB the¬ 
ory with a linear coupling function was considered, and 
the early-time limit of the complete solution of the set of 
field equations was derived; it was found to be identical to 
the solution derived in the context of the pure scalar-GB 
theory. As expected, the higher-curvature - quadratic - 
GB term dominates over the linear Ricci scalar and, in 
conjunction to the scalar field, determines the form of the 
cosmological solution at early times when the curvature 
is strong. That is, in the very early universe it would be 
higher-order curvature terms that would be dominant; 
therefore, we should have a theory involving higher or¬ 
ders of Riemann curvature, and yet one that should be 
ghost-free. This picks out pure Lovelock theory, which is 
a homogeneous polynomial of degree N in Riemann ten¬ 
sor where linear = 1 is Einstein and quadratic N = 2 
is GB. 

Guided by the aforementioned results, in Section 5 we 
proceeded to study the pure scalar-GB theory this time 
with a quadratic coupling function between the scalar 
field and the GB term. Although of a simple form, this 
choice led to a plethora of solutions with interesting char¬ 
acteristics: for a negative coupling parameter, the set of 


equations supported solutions that were either inflation¬ 
ary, de Sitter-type or more involved expanding solutions 
with a de Sitter phase in their past and a natural exit 
mechanism at later times; for a positive coupling func¬ 
tion, the set of equations gave rise to singularity-free so¬ 
lutions with no Big-Bang singularity. All these solutions 
were derived in an analytical, elegant way that allowed 
for the comprehensive study of their properties instead 
of the numerical study that is usually employed in the 
context of similar generalised gravitational theories. 

The case of the general polynomial coupling function 
was briefly considered in Section 6. There, it was demon¬ 
strated that inflationary, de Sitter-type solutions arise 
indeed only for the case of a quadratic coupling func¬ 
tion and for no other. The asymptotic behaviour of the 
scale factor of the universe was derived, and shown that 
for large values of a, the universe adopts a universal be¬ 
haviour independently of the exact form of the polyno¬ 
mial coupling function; for small values of a, on the other 
hand, the form of the scale factor depends strongly on 
the value of the integer n, that determines the power 
of the polynomial function: again, in the strong cur¬ 
vature regime, it is only the case of n = 2 that leads 
to singularity-free solutions whereas all the other choices 
lead to solutions with either an initial or a final singular¬ 
ity. 


Naturally, as the universe expands, its curvature be¬ 
comes smaller and the role of the Ricci scalar will grad¬ 
ually start being of importance again. A toy transitory 
solution was considered in Section 7, that helps to visu¬ 
alise how the transitions between the different eras take 
place: starting from the era where the Ricci scalar can 
be totally ignored, passing to the era where it should be 
taken into account and finally arriving at the one where 
the Ricci scalar dominates thus restoring the traditional 
cosmology. 

Our analysis has singled out the quadratic coupling 
function of the scalar field to the Gauss-Bonnet term as 
one of particular importance: it is only for this choice 
that inflationary solutions as well as singularity-free so¬ 
lutions - two classes of solutions of significant interest in 
cosmology - emerge in the context of the theory. We do 
not consider this to be accidental for the following two 
reasons: on the one hand, the quadratic inflaton potential 
is the only polynomial form that is still compatible with 
the current observational constraints 31|; in our infla¬ 
tionary solutions, the inflaton effective potential is again 
quadratic and, although similar to the one of chaotic in¬ 
flation, it avoids the caveat of the trans-planckian initial 
conditions 2^ . It is worth noting that this attractive in¬ 
flationary model follows from a legitimate string-inspired 
theory rather than being built in an artificial and ad hoc 
way as is usually seen in the literature. On the other 
hand, our quadratic coupling function belongs to a more 
general class of coupling functions that share a number 
of common characteristics with the exact coupling func¬ 
tion of the heterotic string effective theory. It was in the 
context of the latter theory that singularity-free solutions 
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first emerged and the same achievement was performed 
also in the present analysis. We are thus led to conclude 
that a more fundamental, underlying connection may ex¬ 
ist between the quadratic coupling function and the one 
of the ultimate theory, a connection that certainly needs 
to be investigated further. 
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